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1. Overview: Light that is new 

In the mid 1990s the author developed an algorithm that approximates 
the roots of one-variable polynomials in degree six. The method employed 
the dynamical properties of a special map h on CP^. A distinguishing fea- 
ture of this approach to equation-solving is that h is equivariant; it respects 
the symmetries of the Valentiner group V — a collection of projective trans- 
formations isomorphic to the alternating group Aq. Recall that Aq is the 
Galois group of a generic sextic after extracting the square-root of its dis- 
criminant. Being equivariant, h sends a group orbit to a group orbit. In 
algebraic terms, h commutes with each element of V. While the geometric 
behavior of h is elegant, its crucial global dynamical properties were and 
remain unproven. In particular, it was conjectured that a certain Valentiner 
orbit of 72 superattracting periodic points has basins of attraction that form 
a dense subset of CP^. 

Recent results of work on a certain family of complex reflection groups 
[Crass 2005] as well as on an iterative solution to certain seventh-degree 
equations [Crass 2007] led to a substantial dynamical advance: the discovery 
of a new "V-equivariant" g whose dynamics is known to have the desired 
global properties. A reflection group — a class to which the Valentiner action 
belongs — is generated by transformations each of which pointwise fixes a 
reflection hyperplane. [Shephard-Todd, 1954| Indeed, g fits into a pattern 
seen in a number of reflection groups. Namely, it both fixes and is critical 
on the reflection hyperplanes — lines, in this case. This property of being 
critically finite strongly influences the map's dynamics. So much so that 
proofs for its global dynamical behavior are no longer elusive. 

The following discussion 1) examines the geometry and algebra of the 
Valentiner group that is essential for constructing and analyzing g, 2) ex- 
plicitly computes g and establishes its crucial properties, and 3) harnesses 
g's dynamics to a reliable iterative algorithm for solving the sextic. 'Re- 
liable' means that the procedure succeeds in producing a root on a open 
dense set's worth of polynomials and initial conditions. 

Computational results are due to Mathematica. Dynamical plots were 
generated by Julia [McMullen 2011] and Dynamics 2 [Nusse-Yorke 1998] . 
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2. Elements of the Valentiner group: Mirrors, polynomials, 

AND MAPS 

The treatment here is minimal, pointing out only those features relevant 
to the construction of g, the special V-equivariant in degree 31. More thor- 
oughgoing examinations of the Valentiner group appear in jCr ass-Doyle 1997 
and, especially, [Crass 'T999j . 

2.1. Lifting the action. Working with the Valentiner group computation- 
ally requires having a matrix representation. The minimal lift of V to 
occurs in a one-to-three fashion yielding the group V3 where detT = 1 for 
all T G V3. Extending to 

Ve = V3 U-V3 

gives a complex reflection group. (See [Shephard-Todd, 1954| .) 

2.2. Basic geometry. Each of the 45 transposition-pairs {ab){cd) in Aq 
corresponds to a complex reflection on CP^ that belongs to V. Such a 
transformation fixes a CP^ pointwise — a complex mirror — as well as a point 
not on the mirror. For convenience, call these objects "45-lines" and "45- 
points" respectively. 

Among alternating groups, is distinguished by having two systems of 
six versions of ^5 — the six things that Aq permutes under conjugation. The 
Valentiner counterpart Za or 2c (o, c = 1 ... 6) to an A^ subgroup preserves 
a quadratic form Ca or Cg and, hence, the associated conies 

C„ = {C, = 0} or Cc = {Cg = 0}. 

Each conic is, by parametrization, a CP^ thereby endowed with an icosa- 
hedral structure. 

Along with the 45-lines, the icosahedral conies are fundamental to Valen- 
tiner geometry. Another basic piece of the Valentiner structure concerns 
the special V-orbits where 45-lines intersect each other as well as the conies. 
Five of these lines meet at a 36-point, four at a 45-point, and three at a 
60-point associated with either the "barred" or "unbarred" family of conies. 
Special points and certain icosahedral groups have a taxonomic association 
that gives rise to a convenient terminology. A pair of barred or unbarred 
icosahedral conies determines an octahedral structure, namely, 

Oac = St{Ca\JCc)^Si and = St(Cb U Q) ~ ^4 

are octahedral groups by virtue of their preserving the conies 

Cac = {Ca + Cc = ^} and CM = {Cb + Cd = Q}. 

Here, St(X) denotes the stabilizer of X. Table [T] summarizes the special 
structure. 

It turns out that there is an antiholomorphic involution B that exchanges 
the two systems of icosahedral conies. A choice of coordinates yields B{y) = 
y so that B acts as a kind of reflection through an RP^ — the points with real 
y coordinates. Indeed, there are 36 such antilinear maps Bab ™ — one for 
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each barred-unbarred pair — that form a set of generators for the extension 
group 

V = V U SV. 
Of course, each B^b fixes an RP^ mirror Tlab- 

2.3. Generating invariants and equivariants. Under Ve, the forms Ca 
(or Cc) are permuted up to a factor of a cube root of unity. Thus, the 
(or C3) undergo simple permutation; that is, 

for all T e Ve and c G {1, . . . , 6}, o T = C^(„) where r(a) G {1, . . . , 6}. 
Accordingly, 

6 6 
a=l c=l 

is Ve-invariant (a is a simplifying constant). Moreover, F is the Ve invariant 
of least degree and, from it, we can derive two other Ve forms 

$ = det H{F) ^ = det BH{F, $) 

where H(F) is the hessian of F and BH{F, $) is the bordered hessian of F 
and Furthermore, the algebraically independent set {F, generates 
the ring of Ve-invariant polynomials C[y]^''. 

Another form that plays an important role is the jacobian determinant 

X = det J(F,$,^) 

which factors into 45 linear forms that define the 45-lines. However, while 
X is Va-invariant, it enjoys only relative Ve invariance. Specifically, 

XoT = -X for T € Ve - V3 

so that is Ve-invariant and thereby expressible in terms of i*", and ^. 
We can form Ve-equivariant maps by manipulation of the basic invariants: 

^ = VF X V$ (?!) = V-F X V* / = V$ X V* 

where V indicates a formal gradient operator. To obtain a set of generators 
for the V3 equivariants as a module over C[y]^^, include three more maps: 
the identity, 
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where /i64 and k-jo are Ve equivariants of degree 64 and 70 — expressed in 
terms of 'i/'j <P, and / — that vanish on {X = 0}. After division by X, the 
result is a degree-19 map h ("19-map" for short) and 25-map k. 

3. The 31-map: Algebra, geometry, and dynamics 

3.1. Algebraic derivation and geometric properties. Consider the fam- 
ily of 76-maps, 

/76 = (ai + a2 F^^ + ■ ■ ■ + aw^^) i; 

+ (6i • • • + 66 + (ci + • • • + C5 F^) f. 

Utilizing the technique that produced the maps h and k, determine coef- 
ficients so that X divides each coordinate polynomial of fjQ and, thereby, 
obtain a seven-parameter family of 31-maps 

/31 = 

Next, by solving a nonlinear system of equations, force /31 to be critical on 
a 45-line and, thus, on all such lines. The resulting map g has a number of 
special geometric and dynamical properties that make it particularly well- 
suited for residing at the core of a sextic-solving algorithm. First, consider 
g''s algebraic form. For later purposes, having its degree-76 expression is 
convenient: 

X ■g = -162 (67 + 812 ^F"* + 699 ^^F^ - 741 - 150 $^)/ 

- 18 (19 - 1016 - 13457 ^^F^ + 18 ^F^ - 11970 $^F + 13500 
+ (-76 Fi° - 564 $F^ + 5244 ^^F^ - 48672 ^F^ + 66028 ^^F" 

- 591336 ^^-F^ + 232848 ^^F^ - 517752 $2*F + 70200 + 577125 
Since g is critical on Af = {X = 0}, X divides g's critical polynomial 

Cg = det 

where Jg is the Jacobian matrix of a lift ^ of 5 to C^. Moreover, the degree 
of Cg is 90 = 3(31 — 1) so that also divides Cg. That is, for some non-zero 

Cg=l3X^. 

Adapting arguments from [Crass 2005] establishes some basic properties of 
9- 

Theorem 3.1. In the family of V-equivariant 31-maps, g is the unique 
holomorphic member that is doubly- critical on X . 

Proof. Let ^ be a lift of g to and make a linear change of coordinates 

{yi,y2,y3) = r(xi,x2,x3) 

so that L = {2/3 = 0}. Denoting the map in y also by g, 

2/2,0) = (s(yi,y2),i(yi,y2),0). 
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If g{a) = (0,0,0), the homogeneous components s and t have a common 
factor which, contrary to fact, requires the resultant of s and t to vanish. □ 

3.2. Dynamical behavior. The central dynamical object here is the in- 
variant Fatou set J^g where {5^} is a normal family. Its complement is the 
Julia set J'g. 

Since a 45-line L is fixed pointwise by an element of V, a holomorphic 
V-equivariant fixes L setwise. Accordingly, g is strictly critically finite (see 



Fornaess-Sibony 1994||). An intersection of 45-lines is a fixed point that is 



superattracting both along and away from the line. Elsewhere on a 45-line, 
g is critical in the off-line direction so that the line is a superattracting set. 
Following [Crass 2005j and [Ueno 2007] , we can see that g possesses two 
global properties — what 'reliable dynamics' means. The arguments' details 
can be found in these sources. 

First, a computation that expresses g when restricted to a 45-line L: in 
symmetrical coordinates where the four 45-points on L are 

^1/4 



where 

p{z) = - 



{60~^/^(±l±i)}, 

9\l{z) = —-T 
9^) 



- 15(60381703125z3° - 87218015625^2® + 2196104906252^6 

- 169714828125^24 + 21618140625^^2 - 52377193125^20 + 105788081252^® 

- 145978496252^6 - 555500625^^^ - 105496875^^2 _^ 6315475z^° 

- 20343752® - 2205652^ - 57352^ - 46522 - 3) 

qiz) = 2(1153300781252^°- 119174414062522® + 97987851562522^ 

- 2512373203125224 + 1544853515625222 + 31972092187522° 
+ 3560519531252^® - 1249876406252^6 + 2189677443752^ 

+ 105788081252^2 _^ 34918128752^° + 960806252® + 502858752^ 
+ 43379852^ + 11485522 + 5301). 

The remaining intersections of 45-lines — 36, 60, and 60-points and L are 
identified in Figure [TJ 

Now, we consider the content of g^s Fatou set. 

Theorem 3.2. The 36, 45, 60, and 60-points are the only superattracting 
fixed points of g. Moreover, their basins exhaust the Fatou set J-g of g. 

Proof. Since the 45-lines don't blow down to points, their intersections are 
fixed and critical in all directions. Conversely, if a is critical in all directions, 
it belongs to a 45-line L and is a critical point under gi = g\L. But, the 
sixteen points of intersection between L and the remaining 45-lines account 
for all such points. To see this, compute the polynomial that determines the 
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Figure 1. Critical points for g restricted to a 45-line 



critical points oi ql- 

p'{z) q{z) - p{z) q'{z) = a a{zf c{zf d{zf 



where a G C and 



a{z) = ISz"^ + 30^2 - 1 
h{z) = l'az'^ + 1 



c{z) = 15 z*- ( 1 



d{z) = 15/- (^10 + 8 y-ij z^-1. 

Expressed in the z coordinate on L, the 36, 45, 60, and 60-points on L are 
the four simple roots of a, b, c, and d respectively. The multiplicity accounts 
for 

4 (5 + 4 + 3 + 3) = 60 
critical points, exactly the number that a degree-31 map on CP^ has. 
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As for the claim regarding the Fatou set, let F be a component of J-g. A 
nice result in multi- variable complex dynamics states that a strictly critically- 
finite holomorphic map / on CP^ with an irreducible component A in its 
critical set such that f{A) = A has a Fatou set consisting only of basins of 
superattracting periodic points. ( [Fornaess-Sibony 1994| , Theorem 7.8) We 
already have that g satisfies these conditions. Since the only superattracting 
periodic points under g are the intersections of 45-lines, F belongs to a basin 
of one of these fixed points. □ 

Next, we take up a question of the topology of the Fatou set. 

Theorem 3.3. The basins of the superattracting fixed points are dense in 

CP2. 

Proof. Take J & Jg and let U hea, neighborhood of j. By a multi-dimensional 
analogue of Montel's Theorem, U contains critical or pre-critical points, that 
is, 

oo 

u g-\x) n / 0- 

fc=o 

(See [Crass 2005j . Theorem 5.5.) Consequently, some g"^{U) eventually 
meets X and hence a 45-line L. Since g\L is critically finite with fixed 
critical points, its Fatou set Fg\j^ C L consists of basins of the superattract- 
ing points in L. But, these basins are dense — indeed, have full measure — in 
L. Thus, 

r(f/) n / 

so that g"^{U), hence U, contains points in the basin of a superattracting 
fixed point. (Note that the claim pertaining to the Fatou set of g in Theo- 
rem [32] is a corollary to this result.) □ 

Plotting the basins of attraction on a 45-line L reveals sixteen distinct 
basins, one for each intersection of L with a cluster of 45-lines. Four clusters 
of four lines meet L in 36-points, four clusters of three lines occur at 45- 
points, and eight clusters of two lines determine 60 and 60-points. Rendering 
the CP^ as a sphere in Figure [21 each colored region represents a basin for 
one of the superattracting points. The prominent regions indicate immediate 
basins each of which corresponds to one of the points in Figure [1] projected 
onto the sphere. Although the overall symmetry of L is dihedral (-D4), 
the symmetry on the line is Klein-4, since reflection through L acts as the 
identity on L. 

Due to its V symmetry, g preserves each TZab — the RP^ mirror fixed point- 
wise by one of the anti-holomorphic generators of the extended Valentiner 
group V. Figure [3] shows the basins of attraction for eleven superattracting 
fixed points on the RP^. (Here, the plane is projected onto a disk.) In these 
coordinates, (0, 0) is the 36-point to which the RP^ corresponds while five 
36-points and five 45-points are placed symmetrically about the center. In 
fact, the 45-points are on the line at infinity — the boundary circle. Here, the 
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plane's D5 structure is apparent with the five hues of reflective symmetry 
occurring where certain 45-lines meet the real plane in an RP^. Two of 
these RP^s show up as lines of reflective symmetry in Figure [2j 




Figure 2. Basins of attraction on a 45-line. (Produced by Julia.) 



4. A NEW MAP YIELDS A NEW ALGORITHM 

Although the general form taken by the current sextic-solving procedure 
is not new, some of its specific features differ from previous methods with re- 
spect to practical implementation. (Again, an extensive description appears 
in [Crass 1999j .1 

4.1. Sixth-degree resolvent. Recall the six quadratic forms that one 
system of icosahedral subgroups in V3 preserves, up to a factor of a cube 
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Figure 3. Basins of attraction on a 36-RP^. (Produced by 
Dynamics 2.) 



root of unity. By designating the six values 
as roots, the sextic 



k=l k=l 



has coefficients ak{x) that are Ve-invariant rational functions. Hence, each 
ak is expressible in terms of the Valentiner parameters 

$(x) ^{x) 



F{xY " F{xY 
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giving a family of sextic resolvents 



90 ' ' 24300 

(9(30 + iVl5)Vi + 57iVl5 + 100)z^ 
12301875 

(mVi (3(4z + VT5)Vi + 8\/l5 + 42i) - 17i\/l5 - 152)^2 
2214337500 

(3yi(9(-25 + 33iVlE)Vi + 386zv^ + 150) ~ 1944iVT5y2 + 103j\/l5 + 425)-z 

14946778125000 

9Vi(9Vi ((45 + niVT5)Vi - iVTE + 25) - 7iVl5 + 15) - 3i\/l5 - 5 



2421378056250000 

Hence, each is expressible in terms of the Valentiner parameters 



F(x)2 ^ F(x)5 

giving a family of sextic resolvents. After making a simplifying change of 
parameter, the family is given by 

_ -U ^3 

Pviz) 



z(3\/l5- 5i)(36Vi + 1)^ 



■ z 



2421378056250000 
(1296(65 - iVTE) + 144(130 + 21^715) Vi + V2 
14946778125000 
(162(19 + SiVTE) + (-54 + 594i\/l5) Vi - UiVTE - 152) 2 
^ 2214337500 ^ 



(27(35 - 9iV15) Vi + 57z V15 + 100) 3 

12301875 ^ 

(-72^1 - 1H\/T5 + 11) 4 

+ ^ z 

24300 

+ ^i(yT5 + 5i) ^5 + ^6 

The crucial observation is that the 31-map g provides a symmetry-breaking 
tool for solving members of Py. 

4.2. Parametrizing the Valentiner group. Define the map 

y = TxW = ^(x) xwi + $(x) h{x) W2 + F{x) k{x) W3 

that is linear in w = {wi,W2,W3) and degrcc-31 in the parameter x = 
{xi,X2,X2). Parametrize in x a family of Valentiner groups 

where V is the Valentiner group in a selected system of coordinates y. (Recall 
that h and k are V-equivariants of degrees 19 and 25 respectively.) 

Accordingly, an x- parametrized version Fx{w) = F{Ty.w) of the invariant 
F arises. Moreover, each coefficient in Fx{w) of the monomials in Wk is a 
Vx invariant of degree 6-31 and is, thereby, expressible in terms of the 
basic invariants F{x), ^{x), and ^'(x). When normalized to a V-invariant 
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rational function of degree zero, the coefficients are given in terms of the 
parameters V = (Vi, V2): 

= fW 

Similar considerations yield parametrized forms $y(t(;), "^viw), and Xy{w). 

From the ^/-parametrized basic invariants, generating equivariant maps 
tpv, fv follow as does the critically-finite 31-map gy- (For an exactly 
analogous derivation of gv, see [Crass 2007] . To this end, the degree- 76 
expression X ■ g is especially handy.) 



4.3. Finding a 45-line. For a random initial point po in the tf-space, the 
dynamical "output" 

9viPo) P 

is an intersection of 45-lines in the CP'^ under the Valentiner group corre- 
sponding to the parameter choice V. At this point we encounter a complica- 
tion not seen in the original sextic-solving algorithm: the possible values of p 
are not members of a single Valentiner orbit. Furthermore, the combinato- 
rial properties of the various intersections make for an awkward procedure 
that employs the "point datum" p. For instance, a 36-point (at five 45- 
lines) is canonically associated with one conic — hence, one root of Py — in 
each system, whereas a 45-point (where four 45-lines meet) is associated 
with a pair of conies in each system. In order to develop a procedure that 
avoids this pitfall, we can use a dramatic property of the 16-map i/j whereby 
each 45-line collapses to its companion 45-point. 

The idea is to iterate not until a fixed point is located, but only until a 
45-line is found. Recall that each 45-line is superattracting in the normal 
direction so that pk = gyipo) will first feel the influence of a line L and 
subsequently that of a superattracting fixed point on L. Say the trajectory 
Pk has converged "linewise" but not pointwise so that 

Lk = apk-i + Ppk 

is nearly a 45-line. Application of ipy produces a near 45-point 

To achieve high precision we can now compute the trajectory {gv{o,kY \ i = 
1 . . . } to convergence at a 45-point a, a result that provides for two roots of 
Pv. 
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4.4. Selecting a pair of roots. By direct computation, there are Vg- 
invariant and degree-0 rational functions Ak{y) and Bki{y) with the fol- 
lowing properties: 

./ ^ \lk = boTk = d 
^KFacbdj |o A; / 6 and A; / d 

„ , , \ 1 k = b and i = d , , 

BkeiPacbd) = , l<k<i<6. 
10 k or t a 

In the usual way, the functions 

6 

Aviw) = '^Ak{Txw) Sk{x) By{w)= ^ BkiiT^w) SkSi 

k=l l<k<l<G 

have Vx-invariant w coefficients and so, can be F-parametrized. Evaluated 
at a 45-point 

in the t(;-space, 

A = Ay (P45) = sh + Sd and B = Byip^^) = Sb Sd- 
Now, compute the roots Sb and Sd by solving the quadratic equation 

z'^ - Az + B = 0. 

For almost any choice of V and initial point pQ, the dynamics of gy finds 
a 45-line — hence, a 45-point — in the w-spa.ce. Accordingly, a pair of roots to 
Pv is computable from the iterative results. (A Mathematica notebook with 
supporting files that implements the algorithm is available at [Crass web] .) 
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